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1 Introduction 


The field of artificial life, or simply “alife” as it has become known to its 
practitioners, owes its existence and its development to Christopher Lang- 
ton. As Levy's pop history of alife documents [8], Langton picking up the 
cold, dormant trail left by Ulam, Von Neumann, and to a lesser extent Con- ` 
way and the “Game of Life” cult, created an easy to fathom self-reproducing 
cellular automaton and, to underscore its significance, then proceeded to ax- 
iomatize principles for artificial life. Further, under-the auspices of the Santa 
Fe Institute for Nonlinear Studies, Langton organized interdisciplinary con- 
ferences which served both to identify alife investigators and to collect and 
organize their body of work in invaluable, if lengthy, proceedings [6] [7]. 

To over simplify, the rationale for alife runs as follows: If one wants to 
study the origins and development of life, there is one and only one example 
to study — earthbound, carbon based, DNA inspired biological life, known 
as wet life. To find other examples one must look to artificial or synthetic 
life, usually realized in terms of computer simulations through which one 
hopes to distill, isolate, and control the basic tenets of growth, metabolism, 
evolution, behavior, etc. Thus artificial life becomes the antithesis of ar- 
tificial intelligence — the former trying to identify the simplest and most 
elementary components necessary for complex processes while the latter tries 
to mimic, even surpass, complex processes. As one wag has put it, “Why 
should we expect to be able to build a robot with more functionality than a 
human when we can't build a robot with more functionality than an earth- 
worm?” Our purpose is not to survey the scope and history of alife, the, 
reader is referred to the Alife Proceedings for that, so we shall begin to 
narrow our focus by considering the genotype problems alife confronts. 

The point of view adopted by many alife researchers is that since genes 
determine traits, an alife implementation need only specify the traits to 
be studied and then rely on data structures to maintain and modify those 
traits. An obvious way to do this is to view the traits as parameters, and 
an arguably sound way to parameterize traits is via (bit) strings. Fortu- 
nately, or unfortunately as the case may be, this allows the alife researcher 
to use the tools developed by researchers in the field of genetic algorithms, 
a subdiscipline of computer science. 

Genetic algorithms automate algorithm development by evolving and 
evaluating algorithms on the basis of their ability to perform a specific task. 
Again, a survey of genetic algorithms is beyond our scope, but [10] is a com- 
prehensive up to date reference. The point is, that from genetic algorithms 


one can usurp a host of techniques for structuring and evolving bit strings 
that are treated as a genotype of parameters. To furnish an example, we 
shift our attention to computer graphics. 

Two fundamental problems the alife researcher must confront are the 
problem of understanding what is actually taking place during the alife sim- 
ulation, and the problem of interpreting the results of the simulation. How 
does one collect and make sense of data from a simulation that encounters 
billions of “organisms” per run? Scientific visualization and supercomput- 
ers are a must. The most publicized and spectacular alife experiments (for 
example those of Hillis or Jefferson as recounted by Levy [8]) are vivid proof 
of this. The trail we want to follow, blazed by Dawkins [1], combines the 
tasks of genotype representation together with understanding the simulation ` 
by making the organism itself a visual entity! "The details are of interest. 
Dawkins posits a genotype with nine genes controlling embryological traits 
with names like “segment-number” and “segment-distance.” The gene val- 
ues are decimal integers, so the genotype is a decimal string, and the gene 
values drive a recursive algorithm that generates the visual organism or 
biomorph. It is straightforward to observe how perturbing the genotype 
string “evolves” the organism. 

Dawkins? seminal notion of using a genotype string to drive a rendering 
algorithm has been exploited on graphics merits alone. Todd, providing 
technical expertise, and Latham, providing artistic direction, developed their 
Mutator System for allowing an artist to evolve textured and ray-traced 3D 
organisms through a user interface that directly manipulates the gene strings 
[12] [13]. Since there are 24 parameterized genes to control, one thinks of 
the artist as selectively evolving organisms while exploring or journeying 
through 24 dimensional parameter space. Another instantiation of this idea, 
this time for sculpture, is described by McGuire [9]. 

A departure from fixed length, parameterized genotype strings for creat- 
ing visual organisms was developed by Sims [11]. The technique, described 
more fully in the next section, was again implemented in a system whereby | 
users selectively control the evolution of organisms. Thus neither Todd and 
Latham nor Sims ever attempted a true alife simulation using their visual 
entities! We seek to remedy this oversight by describing a series of alife 
experiments with visual organisms designed according to a variant of Sims? 
genotype representation. Before doing so we should mention the one true 
alife simulation whose visual aspects influenced us. 

Larry Yaeger's alife simulation running on a Silicon Graphics Iris Work- 
station is called Poly World. Our description of Poly World is based on notes 


taken during an oral presentation and video demonstration given in the 
Artificial Life Panel Session of SIGGRAPH "92. In PolyWorld the visual or- 
ganisms roam on a bounded two dimensional grid. The organisms “brains” ` 
are small neural nets enabling the organisms to control their external visual 
appearance and to perceive the external world by processing pixmaps. The 
simulation controls for total energy while striving to explore competition 
and self-organization. Genes present are for size, strength, maximum speed, 
mutation rate, number of crossover points in the neural net, lifespan, en- 
ergy to offspring, and ID (a parameter used to enable mimicry). The neural 
net can make decisions about whether the organism should eat, fight, mate, 
move, turn, light (effecting the external appearance of the light sensor panel 
it emits), or focus (gaze at the appearance of others). To see the organisms 
evolve to different species adopting distinct and atypical strategies and be- 
haviors for survival is most impressive. Words do not do justice to the video 
animation sequences. 


2 Graphical Alife 


In this section we describe our implementation of Sims’ mutating expressions 
[11]. The backbone of our implementation was first developed for a computer 
art application [3]. The important differences to note in our implementation 
are: (1) Where Sims uses neighborhood image processing functions (e.g.,“ 
warping, blurring, dissolving, filtering, and noise generating functions) we 
use point processing functions which are always normalized so that inputs 
and outputs lie in the unit interval and (2) When visualized, our organisms 
are rendered relative to a cmax valued color map as opposed to an absolute 
rgb scale. Unfortunately, the sequence of experiments described in the next 
section did not progress to the point where environment and organisms were 
rendered in color. 


Let V be a set of algebraically independent variables, let J be the unit 
interval, and for i = 1,...,n fix an integer d; > 0, a function f; : 1% — I, 


and let F = {f;}. 


DEFINITION. An ezpression is defined recursively as a node (a,b, g) 
where a,b € I, and either g € Y or g € F. When g E F, the arguments of g 
are again nodes. 


We will use V = {c, u, v, t}, where c represents a constant, u and v param- 
eterize the unit square, and t is a time variable. We only require mappings 


f: defined on I, 1?, and eventually I? (viz., d; = 1, 2, and eventually 3). 
Thus, until further notice F consists of functions from the unit interval to 
itself with one argument named 


Ínsin, ncos, nexp, nlog, nabs, nsqt, nsqr, ncub, nneg) 


together with functions from the unit square to the unit interval with two 
arguments named 


{nadd, nmul, nmod, nmin, nmax, npwr, nand, nvee, ncir}. 


The precise definitions for these functions are not relevant, but the mnemon- 
ics will give clues about their construction. 


EXAMPLE. Omitting outermost parentheses for expressions, suppressing 
commas between arguments, and remembering that nsin, ncos, and nsqt 
are suitably normalized sine, cosine, and square root functions in one vari- 
able, a legal expression which up to coefficients calculates the minimum of 
the sine of u and the cosine of the square root of v is 


0.05 0.38 nmin(0.54 0.09 nsin(0.42 0.52 u) 0.03 0.69 
ncos(0.22 0.40 nsqt(0.04 0.63 v))). 


DEFINITION. An expression is evaluated at (to, Uo, Vy) using the recursive 
evaluation rule € given by 


čla,b,g) >axélg) +b (mod 1) 


where if g € Y, 


0 ifg=c 
£g) = uo ifg=u 
vy ifg=v 
and if g = f(ni,...,na), where f € F and n,..., ng are the argument nodes 


of f, 
E(g) = f(E(m),..-,E(na)). 


To summarize, our organisms are computational entities that can also 
be visualized. We next describe their evolutionary and mating capabilities. 
A “creation” event provides a pseudorandom expression, used as a starting 
genotype, that is very primitive, perhaps something as mundane as 


0.85 0.38 ncir(0.80 0.35 v nneg(0.42 0.52 u)). 


This primitive genotype is “evolved” to provide the complex genotype for a 
species. This process is accomplished by successively feeding the genotype 
to an algorithm that complexifies it by evolving nodes of type V into those 
of type F on a probabilistic basis, and similarly evolving some one argu- 
ment nodes into two argument nodes. The probablistic genotype evolution 
algorithm has many decision table parameters and is “intelligent” enough , 
to try and preserve the history of the genotype by promoting evolution near 
the leaves rather than near the root of the expression. A typical species 
genotype that might arise follows. . 


0.87 0.56 nmul(0.80 0.63 nsin(0.85 0.95 nsqt(0.80 
0.62 u)) 0.87 0.75 nexp(0.80 0.17 nmin(0.80 0.40 v 
0.80 0.47 u))) 


All organisms of the species will use this genotype, meaning their expressions 
will have this identical tree structure. 

The initial population of the species consists of NORGANISMS, a third of 
which have all their nodes randomly selected from the V and J sets while 
still remaining faithful to the the genotype; a third of which have a randomly 
selected proportion of their nodes determined in this fashion; and a third 
third of which have all arguments faithful to the species genotype but with 
the a and b coefficients at each node severely perturbed. Our alife simulation 
follows the development of this population, or colony, through a prescribed 
number of time cycles, say NCYCLES. 

The organisms can reproduce sexually or asexually. The sexual reproduc- 
tion, or mating, algorithm for two organisms creates an offspring organism 
faithful to the species genotype node by node with reference to a parameter “ 
that biases node selection in favor of a dominant organism. Asexual repro- 
duction merely perturbs the node's a and b coefficients as described above. 
All organisms are subject to mutation which can alter the árguments of 
one or more nodes, and all organisms undergo “genetic drift” which is im- 
plemented by a very mild perturbation of each node's a and b coefficients. 
There are many more aspects of the simulation to discuss, but they will be 
dealt with as they arise in the context of the sequence of experiments that 
was actually performed. 


3 The Experiments 


Experiment #1. 


In general, our organisms are camouflage organisms whose survival de- 
pends on their ability to blend in with their surroundings. For this first 
experiment, we challenged them to breed to match a simulated background. 
Since the organisms perform computations on the unit square via the u,v 
parameters, we assigned random values r, , to the following nine “sites” of 
the unit sguare: 


[(u;, 05): u; = 1/3, v; = 7/3,0 < i,j < 2). 


We used the sum of the squares of the errors to these target values as the 


fitness criterion: 
C(g) = X (E(g) (use) — Ti)“: 
ij 


é 


The population was held steady at 100 organisms. After each time cycle, 
the best 20 organisms were retained for breeding. The most fit organism 
was mated with the other 19, and the second most fit organism was mated 
with the other 18. The population was then restored to 100 organisms by 
allowing all the breeding stock to contribute aserually produced offspring. 
Note that the time variable t was not used in our early experiments. 

We have mentioned that every organism is subject to genetic drift and 
mutation. When an organism is created a random number is compared to 
the mutation rate, mrate, to decide whether or not a randomly selected node 
should be altered. The mutation rate begins at 0.20 but is adjusted during 
the course of the simulation so that it may drop as low as 0.02 or soar as 
high as 0.40. At the completion of each time cycle, the rate changes by 0.01 
according as the most fit individual present — the one with the greatest 
influence on the next generation — is better at camouflage than the most fit 
individual of the previous generation. Our typical run was set for 100 time 
cycles. 

We observed that by the end of a run the most fit organism achieved a 
fitness of about 0.05 and the mutation rate hovered near 0.10. Surprisingly, 
almost always, superior organisms were computing using only constants and 
one of the u and v variables! Examining sample sets of target values sug- 
gested an explanation for what was happening. The organisms had taken 
advantage of the fact that many sites had similar target values. With limited 
computational resources an organism was rewarded for focusing its atten- 
tion to the variable most responsible for target variation, and the sum of 
squares fitness favored a genetic attempt at one variable approximation to 
cover this spread. We viewed the breeding stock at the end of each run 


and observed that while node structures were nearly identical with respect 
to arguments, indicating that the breeding stock was inbred, the modest 
genetic drift (a mazimum of 0.05 per coefficient) was severely effecting the 
fitness. Typically the least fit member of the breeding stock had fitness of 
about 0.25. We expected tighter clustering of the fitness values. 

Our simulation does not allow node coefficients to “wrap,” so an a coef- 
ficient drifting toward 1.00 cannot wrap to 0.00 and an a coefficient drifting 
towards 0.00 cannot wrap to 1.00. Node coefficients, especially dominant 
coefficients close to the root of the expression, had an annoying tendency 
to drift to extremes — 0.00 for a and 1.00 for b, or conversely — indicating 
that successful strategies often resulted from ignoring variables or ignoring 
constants. Again the distribution of target values was critical. 


Experiment #2. 

To promote greater diversity in the population and to try and counter- 
act how quickly the simulation was reaching steady state, we increased the 
complexity of the genotype by feeding it through the genotype evolution 
algorithm 20 times as opposed to 10. We increased the initial mutation rate 
from 0.20 to 0.33, we eliminated the dominance factor in sexual mating, . 
and we permitted more dramatic mutations by altering two node arguments 
whenever mutations were called for. But the key change was to view organ- 
ism camouflage sites not as fixed points but as areas which would respond 
to point stimuli. Thus an organism was envisioned as trying to camouflage 
itself by matching as best it could random values assigned to each of the 
nine 1/3 x 1/3 uniformly sized patches or sites of which it was made up. To 
evaluate fitness, a random point was selected in each patch, and its compu- 
tational value was compared to the target value for that patch. The random 
points clustered around the center of the site so the organism had a better 
chance to make graceful transitions from site to site. 

We indeed found that organisms had a much more difficult time with 
the camouflage task. It was rare to see the mutation rate drop below 0.25 
during the course of a run. The variation in target values coupled with 
being kept off balance with respect to site queries meant that in every time 
cycle the majority of the colony was poorly camouflaged, while a few lucky 
organisms were well camouflaged. We suspected that the latter was due to 
being queried at the “right” set of site points. 


Experiment #3. 
At this point, we implemented the truly unique feature supported by our, 


alife design, an environment for the organisms to live in that could also mu- 
tate and evolve! The environment was also created using the genotype evolu- 
tion algorithm. This meant our simulation bore some semblance to a “Gaia 
model” or perhaps a colony parasitic on an individual of another species. 
Hypothesizing that the environment would be perceived as relatively stable 
in comparison with our short lived camouflage species, we only allowed the 
environment to alter itself by genetic drift using the miniscule drift limit 
setting of 0.01. The environment was treated as a 10 X 10 organism. (Recall 
that camouflage organisms are treated as 1 x 1 organisms.) We reduced the 
colony size to 25 and set the mutation rate back to 0.20. The organisms 
now had 25 patches or sites. Randomly selected points within perimeter 
sites were used to calculate camouflage values (after all, why should one be 
able to see the interiors of these two dimensional organisms?), and these 
were now compared to the environment organism's values computed at the 
same points. Obviously, scaling and translation were necessary to “align” 
organism and environment. The x and y coordinates used to determine the 
organism's location were at the organism's discretion because interior sites 
were reserved for their computation. All reproduction was by “cloning” so 
one would expect a newly cloned organism to be located close to its parent 
because its position calculation apparatus would be inherited. After each 
time cycle the best five camouflaged organisms were cloned twice, and the 
worst five were not cloned at all. Henceforth in our experiments the muta- 
tion rate fluctuated at each time cycle according to the average fitness of 
the colony. 

This first experiment using an environment organism was primarily a 
“proof of concept.” The colony camouflage rating and mutation rate quickly 
stabilized, and examining the top ten individuals after each run revealed 
that organisms clustered at the “boundary” of the environment or along the 
diagonal. 


Experiment #4. 

Yet another “proof of concept” experiment. Now we allowed our or- 
ganisms to have a lifespan. We increased the population size to fifty, and 
permitted an organism to live to a maximum age of ten. This was not at 
easy as it sounds. At the end of each time cycle, juveniles, those under age 
two, were “protected,” the least fit adults perished, and the breeding popu- 
lation — the most fit adults aged two through eight — were cloned to bring 
the population back to full strength. The initial population had random age 
assignments. 
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We were pleased to observe that the average age of the colony quickly 
reached steady state. The organisms still clustered, staking out positions on 
the boundaries, on horizontal or vertical strips, and on the diagonal. 


Experiment #5. 

This experiment relaxed the camouflage calculation by querying only 
the four corner sites together with the middle sites at each edge. In future 
experiments we planned to use the remaining perimeter sites for communi- 
cation with the external world. But the primary purpose of this experiment 
was to initiate some movement on the part of the colony by using a random 
point within the designated location sites to (re)calculate the organism's x 
and y coordinates at each time cycle. We increased the size of the colony 
to 70 and the number of time cycles to 75, but the colony remained very 
static and clustered at the diagonal. The fact that organisms clumpeded at 
all was viewed positively — similar to discovering a well camoúflaged patch 
of mushrooms or nest of insects. 


Experiment #6. 

Having assigned random ages to the initial population, we now also ran- 
domly assigned initial locations to the population, and we forced the organ- 
isms to “explore” their environment by using their location sites to compute 
x and y displacements. For the time being we left the decision of whether 
displacements should be positive or negative to chance. The distance per 
time cycle an organism could travel could not exceed its “length.” Newly 
cloned organisms took their place at the side of their parent. Without a 
graphics display it was difficult to understand the simulation, even with as 
few as 20 organisms and 10 time cycles (one lifespan). 


Experiment #7. 

We hooked our alife simulation into the SRGP graphics package [2] al- 
lowing us to view the position of each organism — actually the lower left 
hand corner of each organism superimposed on a display scaled so that the 
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environment was 300 x 300 pixels — after each time cycle. Moreover, mu- * 


tant organisms were specially tagged. Since an organism might traverse its 
environment in a single lifespan, we made the environment toroidal allow- 
ing its offspring to duplicate its migratory path. A computation site was 
designated for the organism to determine a compass direction to use in con- 
junction with its displacements. Here we restricted the directions to the 
four cardinal points of the compass. Organisms were allowed to incorporate 
the time variable into their genotype, but the time variable was only for 
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camouflage organisms not for the “timeless” environment organism. The 
organism's age served as the computational input for time. Also, organ- 
isms were permitted to use “logic” via simple-if, nsif, or if-then-else, nsel, 
nodes. The conditional component of the expression evaluated to TRUE when 
the computational value was greater than one-half, and to FALSE otherwise 
(a kind of fuzzy logic). For simple-if, FALSE meant the organism computa- 
tion yielded 0.0, but in all other cases the natural computational value was 
returned. 

As expected the display revealed that there were zones within the en- 
vironment where the organisms could flourish. But even with movement 
slowed so organisms could move at most half a body length per time cycle 
motion patterns were confusing. 


Experiment #8. 

We slowed the movement so that even under maximum speed in a con- 
sistent direction it would take an organism 50 cycles to traverse its environ- 
ment. We also allowed the organisms to travel diagonally, thus organisms 
could now choose from eight compass directions. We increased the lifespan 
of an organism to 20 cycles, and we re-populated the colony at the end of 
each cycle by mating randomly selected pairs of the breeding population. 
Offspring location was set so that a newborn would accompany one of its 
parents. In this model every member of the breeding population disperses 
its genotype, and some collect this reproductive material for sexual union 
and birth. This is reminiscent of fertilization in tide pools or coral reefs. 
Graphics were enhanced by using distinct markers for juveniles, mutants, 
and adults. 

The simulation was allowed to continue for up to 900 cycles during some 
runs. There was considerable ebb and flow observed in camouflage average 
and in the movement of our colony of 50 organisms. 


Experiment #9. 

We wanted to make the colony more assertive than its environment; 
aggressive learners and evolvers if you will. Thus we designed functions 
supporting logical computation, denoted bsif, bleg, and bgeg for boolean 
simple-if, less-than-or-equal, and greater-than-or-equal. (Warning: The 
nand function mentioned earlier is calculated on the basis of the bit pattern 
of its two arguments). These new boolean two argument functions returned 


d 


fuzzy truth values: for simple-if, a fuzzy FALSE otherwise the computational , 


result; and a fuzzy TRUE or FALSE as necessary for the inequalities. The en- 
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vironment organism was prevented from using the time parameter or these 
fuzzy boolean nodes, but colony organisms cóuld enjoy both. All organ- 
isms could use the boolean if-then-else, now renamed bsel. We also forced 
the organisms to rely more heavily on their computational machinery by 
enforcing a lower bound of 0.80 on the a coefficient at each node. Finally, 
we enlarged the environment to a 100 x 100 world for our 1 X 1 organisms. 
Conceivably, the environment could still be crossed in fifty time cycles. 

We have avoided cluttering experiment descriptions with computer print- 
outs. But appended is a sample graphics display from the end of a 200 time 
cycles run for a colony of 50 organisms. It is accompanied by printout that 
gives the static genotype for the environment, the initial template genotype 
for the colony, and after every tenth time cycle the average age of the colony 
(as a percentage of maximum lifespan), the average camouflage value, and 
the current mutation rate. The printout summarizes this data for every in- 
dividual in the colony at the end of the run, and then provides the genotype 
of every fifth organism in the colony as ranked by camouflage value. 


4 Open Questions, Future Work 


Our alife simulation experiments have raised almost as many questions as 
they have answered. We should first remark that though the genotype of our 
organisms bears a resemblance to that of Koza [4] [5], the inspiration clearly 
comes from Sims [11]. For Koza, the tree structures parse to pseudocode 


a 


and there is no concept of a species genotype, hence all tree structures ` 


are templates for the population. More to the point, Koza’s organisms can 
evolve to a “perfect form” by, say, successfully running a maze or reproducing 
a logic table. We are unable to define perfect organisms. Our use of a 
template genotype seems natural but raises a vexing question. 


Question. How does a species complexify its genotype? Specifically, how 
could our simulation increase the size and complexity of the genotype “on the 
fly” and still avail itself of the node-by-node method for sexual reproduction 
that seems most natural? 


The reason it is important to try and provide answers to these questions 
is because only by complexifying the genotype will a species be able to build 
and pass on a “genetic memory.” It may be seen as a shortcoming of our 
model that our organisms do not “learn” during their lifespan. For now we 
view the organisms as too primitive to act on external stimuli and record 
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information about such actions. We have considered ways to incorporate 
such notions, but to speculate further on them at this time seems wholly 
premature. , 

Clearly, the graphics for our simulation need improvement. First, the 
environment should be visualized so we can “see” what the organisms must 
contend with. Second, we would like to be able to see if the organisms are 
exploiting the fact that exterior sites are for camouflage and interior sites 
are discretionary. Organism and environment visualization would imme- 
diately give clues about the wisdom of using an environment to organism 
length ratio of 10 versus 100. Considerable compute power seems called 
for in the future. Perhaps the simulation should even be performed on a 
supercomputer. There is also a technical visualization question to consider. 


(Juestion. Our organisms are now permitted to use fuzzy logic. Can our 
visualization algorithm cope with this? 


We need to devote more effort to understanding what our organisms 
are computing. Certainly we need tools to visualize their genotypes in tree 
form. This might make their computation strategies apparent or reveal how 
much or little computational control they exhibit. This would also provide 
indications about how the parameters in V are being used, and suggest what 
a reasonable genotype size should be for supporting 25 computational sites. 

Future work should include a series of comparison experiments to explore 
how useful fuzzy logic functions are, and how utilitarian the core set of 
functions F is. It would make sense to try several different species genotypes 
on the same environment in combination with several F designs. We should 
even experiment with other camouflage measures! 


(Juestion. For evolution or alife purposes, how should one best measure 
how well n fitness values e,,...,€, estimate or approximate n target values 


A a, 


Several very hard alife problems are still to be addressed. How large 
should the population be? How long should the simulation run? How do 
we let decisions about reproduction emerge? Several standard models base 
these types of decisions on environmental resource measures and individual 
fitness measures referred to as individual health or energy. We have not 
considered trying to design food consumption and food supply mechanisms 
because we believe camouflage is a different paradigm, one supported by a 
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resource laden environment. In fact, it might make more sense for us to 
design an ecosystem, perhaps consisting of a large, slowly dispersing cam- 
ouflage species and a small, fast moving predator species (think chitons and 
starfish). 

Such a scenario leads us to conclude by re-emphasizing the strengths of 
our model to help justify why we think it has future alife potential: Envi- 
ronment and life forms are governed by a consistent “biology” — mutating 
expressions — that can adapt, mutate, evolve and be set to any task that we 
are willing to formulate as purely computational. — 
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x — adult 
o — mutant 
u - juvenile 
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Al. Organism Distribution 


environs 
> 0.80 0.02 ncir(0.80 0.97 nsar(0.80 0.77 u) 0.80 0.96 
nvee(0.80 0.25 npwr(0.80 0.29 u 0.80 0.71 u) 0.80 0.63 
nneg(0.80 0.84 v))) 


genotype 
> 0.80 0.84 bsel(0.80 0.06 nmin(0.80 0.70 nexp (0.80 

0.35 nadd(0.80 0.62 t 0.93 0.61 nneg(0.80 0.62 t))) 0.89 

0.43 nsgt (0.80 0.12 u)) 0.80 0.01 bsel(0.80 0.67 nvee( 

0.80 0.99 v 0.80 0.72 v) 0.80 0.42 ncir(0.80 0.91 bsel( 

0.80 0.79 u 0.83 0.48 nsgr(0.80 0.30 c) 0.80 0.37 c) 

0.80 0.54 bsel (0.80 0.18 t 0.80 0.67 t 0.80 0.33 t)) 

0.80 0.50 t) 0.80 0.78 nand(0.80 0.90 v 0.80 0.34 v)) 
cycle = 0 tavg = 0.5210 cavg = 1.0849 mrate = 0.19 
cycle = 10 tavg = 0.3400 cavg = 0.5277 mrate = 0.15 
cycle = 20 tavg = 0.3880 cavg = 0.4366 mrate = 0.11 
cycle = 30 tavg = 0.3710 cavg = 0.5516 mrate = 0.13 
cycle = 40 tavg = 0.3970 cavg = 0.4708 mrate = 0.13 
cycle = 50 tavg = 0.3790 cavg = 0.4728 mrate = 0.15 
cycle = 60 tavg = 0.3620 cavg = 0.3439 mrate = 0.09 
cycle = 70 tavg = 0.4040 cavg = 0.5391 mrate = 0.13 
cycle = 80 tavg = 0.4000 cavg = 0.6935 mrate = 0.15 
cycle = 90 tavg = 0.3810 cavg = 0.6123 ° mrate = 0.17 
cycle = 100 tavg = 0.4130 cavg = 0.4337 mrate = 0.13 
cycle = 110 tavg = 0.3960 cavg = 0.4039 mrate = 0.13 
cycle = 120 tavg = 0.3760 cavg = 0.3974 mrate = 0.13 
cycle = 130 tavg = 0.3730 cavg = 0.4217 mrate = 0.13 
cycle = 140 tavg = 0.3830 cavg = 0.2961 mrate = 0.13 
cycle = 150 tavg = 0.3610 cavg = 0.5042 mrate = 0.15 
cycle = 160 tavg = 0.4040 cavg = 0.5060 mrate = 0.17 
cycle = 170 tavg = 0.3160 cavg = 0.7278 mrate = 0.21 
cycle = 180 tavg = 0.4040 cavg = 0.7435 mrate = 0.21 
cycle = 190 tavg = 0.3970 cavg = 0.6708 mrate = 0.21 

A2. Environment and Species Genotype; Simulation Averages 
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Individual Organism Data 
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0.89 bsel (0.80 0.91 ncir(0.80 0.98 nexp(0.95 
nmax(0.82 0.91 t 0.82 0.38 nneg(0.87 0.65 v))) 0.86 
nneg(0.80 0.33 u)) 0.80 0.08 bsel(0.80 0.20 nvee( 
0.97 t 0.90 0.62 v) 0.87 0.59 ncir(0.90 0.15 bsel( 
0.86 u 0.80 0.48 nsgr(0.84 0.63 v) 0.84 0.22 v) 
0.38 bsel(0.80 0.82 t 0.82 0.36 t 0.80 0.07 u)) 
0.76 u) 0.95 0.34 nmod(0.80 0.18 t 0.91 0.82 u)) 


0.96 bsel(0.90 0.12 ncir(0.80 0.81 nexp(0.82 
nmax(0.87 0.72 u 0.86 0.72 nneg(0.86 0.97 t))) 0.80 
nneg(0.90 0.51 u)) 0.80 0.05 bsel(0.92 0.67 nvee( 
0.82 v 0.84 0.79 v) 0.80 0.60 ncir(0.88 0.77 bsel( 
0.56 t 0.92 0.53 nsgr(0.80 0.49 u) 0.80 0.23 c) 
0.80 bsel(0.91 0.99 t 0.80 0.45 c 0.80 0.37 u)) 
0.19 u) 0.83 0.22 nand(0.81 0.80 v 0.80 0.99 u)) 


0.80 bsel (0.86 0.22 ncir(0.83 0.60 nexp (0.87 
nmax(0.85 0.87 u 0.85 0.72 nsgr(0.92 0.66 u))) 0.80 
nsin(0.87 0.89 u)) 0.80 0.16 bsel(0.82 0.72 nvee( 
0.95 v 0.85 0.43 v) 0.94 0.41 ncir(0.96 0.88 bsel( 
0.65 u 0.80 0.23 nsgr(0.80 0.63 u) 0.80 0.30 v) 
0.04 bsel(0.80 0.86 t 0.81 0.78 t 0.80 0.74 v)) 
0.31 u) 0.89 0.10 nmul(0.91 0.90 v -0.87 0.45 v)) 


0.14 bsel (0.88 0.25 nmin(0.97 0.86 nexp (0.80 
nmax(0.86 0.58 t 0.86 0.40 nneg(0.97 0.39 t))) 0.87 
nneg(0.89 0.38 u)) 0.83 0.06 bsel(0.92 0.43 nvee( 
0.70 t 0.96 0.90 v) 0.90 0.53 ncir(0.80 0.17 bsel( 
0.69 u 0.90 0.51 nsar(0.80 0.31 u) 0.80 0.47 v) 
0.29 bsel(0.86 0.92 t 0.89 0.27 u 0.89 0.38 t)) 
0.37 v) 0.84 0.43 nmod(0.83 0.29 t 0.89 0.99 u)) 


0.91 bsel(0.98 0.95 nand(0.80 0.96 nexp(0.85 
nmax(0.80 0.92 t 0.86 0.40 nneg(0.80 0.42 t))) 0.87 
nneg(0.80 0.28 u)) 0.82 0.09 bsel(0.97 0.21 nvee, 
0.90 t 0.80 0.97 v) 0.84 0.50 ncir(0.82 0.69 bsel( 
0.61 u 0.88 0.10 nsgr(0.80 0.69 v) 0.88 0.26 v) 
0.70 bsel(0.98 0.24 t 0.80 0.31 t 0.88 0.28 u)) 
0.47 u) 0.80 0.24 nmul(0.80 0.89 t 0.86 0.75 u)) 


0.97 bsel(0.80 0.10 ncir(0.82 0.56 nexp(0.80 
nmax(0.90 0.48 u 0.81 0.23 nneg(0.80 0.33 t))) 0.80 
nneg(0.80 0.26 u)) 0.86 0.17 bsel(0.83 0.77 bleg( 
0.71 v 0.94 0.54 c) 0.96 0.46 ncir(0.80 0.04 bsel( 
0.58 u 0.83 0.24 ncub(0.83 0.32 u) 0.90 0.20 v) 
0.36 bse1(0:80 0.02 t 0,80 0.53 t 0.84 0.79 v)) 
0.66 u) 0.80 0.40 nadd(0.80 0.06 v 0.80 0.76 u)) 


0.87 bsel(0.93 0.05 ncir(0.80 0.53 nexp(0.80 
nmax(0.90 0.46 u 0.93 0.40 ncub(0.99 0.81 u))) 0.92 
nneg(0.80 0.13 u)) 0.80 0.11 bsel(0.80 0.36 nvee( 
0.52 t 0.96 0.90 v) 0.80 0.46 ncir(0.80 0.79 bsel( 
0.77 u 0.80 0.38 nsgr(0.80 0.50 t) 0.86 0.08 v) 
0.25 bsel(0.80 0.73 t 0.80 0.82 t 0.90 0.81 v)) 
0.55 t) 0.80 0.19 nand(0.88 0.02 v 0.80 0.18 v)) 


A4. Selected Individual Genotypes 


.90 0.76 bsel(0.90 0.35 ncir(0.80 0.59 nexp(0.85 

.19 nmax(0.92 0.38 t 0.80 0.49 nneg(0.82 0.42 t))) 0.90 
.64 nneg(0.82 0.88 u)) 0.89 0.09 bsel(0.80 0.70 bleq( 
0.04 t 0.80 0.32 v) 0.80 0.65 ncir(0.87 0.97 bsel ( 
.80 0.66 u 0.89 0.31 nsgr(0.80 0.66 v) 0.87 0.05 v) 

.83 0.81 bsel(0.80 0.16 t 0.86 0.29 u 0.80 0.22 u)) 

.96 0.67 u) 0.84 0.35 nand(0.96 0.09 t 0.80 0.81 u)) 
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.89 0.99 bsel(0.86 0.12 ncir(0.96 0.62 nexp(0.80 

.33 nmax(0.92 0.52 t 0.80 0.98 nneg(0.80 0.05 t))) 0.80 

.39 nneg(0.96 0.88 u)) 0.96 0.11 bsel(0.80 0.80 nvee( 

0.61 v 0.80 0.61 v) 0.90 0.16 ncir(0.80 0.81 bsel( 
0.66 t 0.87 0.24 nlog(0.95 0.65 v) 0.91 0.25 c) 

.83 0.93 bsel(0.81 0.83 t 0.98 0.25 t 0.81 0.42 t)) 
0.55 u) 0.95 0.45 nand(0.84 0.10 v C.93 0.93 u)) 
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.85 0.04 bsel (0.80 0.06 nmin(0.99 0.74 nexp(0.94 

„44 nmax(0.80 0.17 u 0.80 0.89 nsgr(0.89 0.33 t))) 0.90 
.55 nneg(0.90 0.93 u)) 0.80 0.02 bsel(0.80 0.78 nvee ( 

4 0.94 t 0.98 0.82 v) 0.90 0.39 ncir(0.99 0.86 bsel ( 
.94 0.84 t 0.80 0.39 nsgr(0.89 0.71 u) 0.88 0.31 c) 

.95 0.56 bsel (0.80 0.93 t 0.80 0.82 t 0.80 0.45 t)) 

.90 0.81 u) 0.80 0.50 nand(0.84 0.86 t 0.85 0.74 u)) 
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mutations : 173 
births : 995 


A5. Genotypes continued 


